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Bound states in the continuum in open acoustic resonators 1 

We consider bound states in the continuum (BSC) or embedded trapped modes in two- 
and three-dimensional acoustic axisymmetric duct-cavity structures. We demonstrate 
numerically that under variation of the length of the cavity multiple BSCs occur due to 
the Friedrich-Wintgen two-mode full destructive interference mechanism. The BSCs are 
detected by tracing the resonant widths to the points of the collapse of Fano resonances 
where one of the two resonant modes acquires infinite life-time. It is shown that the 
approach of the acoustic coupled mode theory cast in the truncated form of a two-mode 
approximation allows us to analytically predict the BSC frequencies and shape functions 
to a good accuracy in both two and three dimensions. 
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1. Introduction 

Bound states in the continuum (BSCs) also known as embedded trapped modes are lo¬ 
calized solutions which correspond to discrete eigenvalues coexisting with extended modes 
of a continuous spectrum in resonator-waveguide configurations. The existe n ce of tr apped 
solutions residing in the continuum was first reported bv lvon Neumann fc Wigner ( 192911 ) 


at the dawn of quantum mechanics. To the best of our knowl edge , the term bound 
state (embedded) in the continuum was introduced by iFondal ( 196 1) in the context 


of resonance reactions in the presence continuous channels. Since then bound state in 
the continuum has been u niversally used to designate a BSC in quantum m echanics 
Stillinger fc Herricfl ( 19751 ). In the field of fluid mechanics, IParkerl ( 19661 196 7) is cred¬ 


ited to be the first to encounter resonances of pure acoustic nature in air flow over a 


uide structures ranging from quantum wires Shahbazvan & Rnikh (19941): Kim et al. 

(199fJ); Olendski & Mikhailovska (2002): Sadreev et al. (2006): Cattao; 

an & Lott 

(2007), 

acoustic waveguides Linton & Mclverl (20071): Duan et al. (2007): He 

m & Koch 

(2008); 

Hein et al. (2012), and photonic crystals Shioman & Venakides (2005!); 

Bulgakov fc Sadreev 

(20o 3); Marinica et al. (20081): [Yang et al. (2014). The BSCs are of 

immense interest, 


specifically, in optics thanks to experimental opportunity to confine light in op tica l micro¬ 


cavities despite that outgoin g waves are allowed in the surrounding medium IP lot nik et ali _ 

( 2011 ): ICorrielli ef~al\ ( 20131 ): lWeimann et all ( 20131 ): lchia Wei Hsu et all ( 20131 ) : lz~hang fc Zhand 


( 2015 ). At the same time, in aerodynamics, trapped and nearly trapped modes are known 
to cau s e severe vibrations and noi s e problems in gas and steam pipelines ljungowski et al. 


( 1989 ): Ziada fc Buhlmannl ( 1992f ): lKriesels et al. ( 1995h : Tl 

(2q3- 


Tonon et al. 


The most trivial mechanism of BSC formation is due to the symmetry of the structure 
i.e. the BSC and the continuous-spectrum modes have incompatible symmetries. Typ¬ 
ically, when symmetric and antisymmetric problems are separated the first symmetric 
and antisymmetric propagating modes have different cut-off frequencies. Thus, the dif¬ 
ference between the cut-off frequencies provides a window for an antisymmetric BSC to 
reside in. Su c h B SCs are omnipresent in wave-related set-ups in cluding qu antum wires 


dl989h. fluid in 


a wave 


tank [Evans fc Lint on (Il99ll). acoustic [Evans et all 


: ISugimoto fc Imahoril ( 2005f ) and elastic [Maksimov fc Sadreevl ( 20061 ) waveguides. 


f Email address for correspondence: mdn@tnp.krasn.ru 
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By breaking the axial symmetry of the structure the symmetry protected BSCs be¬ 
come quasi-trapped le a ky mod e s coup l ed to the propagating mc> d es of the waveguides 


I Aslanyan et all (|2000l) . iNockell (|l9921) : lLadron de Guevara et al\ ( 2003f ) demonstrated 


that if the symmetry is broken under variation of a control parameter the symmetry 
protected BSCs reveal themselves in form of ghost Fano resonances. Conversely, as the 
symmetry is recovered the Fano resonance collapses with one of the two resonant modes 
acquiring an infinite live-time. The collapse of Fano resonance is a signature of full de¬ 
structive interfere nce of two degene rate eigenmodes of the resonator on the waveguide- 
resonator interface iKim et all (1 9991 ). In the consequence of this a trapped mode is formed 
as a localized solution decoupled from the extended modes of the waveguides. Impor¬ 
tantly, the full destructive interference of two degenerate eigenmodes_ leaking into waveg¬ 
uides represents a generic mechanism of BSC formation Friedrich fe Wintgen (Il985 ) _ 

whose implementations go far beyond the above symmetry arguments by Ladron de Guevara et al 
( 2QQ3h . In particular, this mechanism allows the formation of B SCs with th e same sym¬ 
metry about the duct axis as the coexisting scattering solution [Hein et al. | (120121) . The 
BSCs which can be attributed t o a fu ll destr uctive interference are mul tiply reported in 
quantum mechan ics Kim et al. (1999), optics Bulgakov fc Sadreevl ( 2014fh and acoustics 
Hein all dgOlcL l2012h. The same interfere nce mechanism was recently employed by 


Lepetit et all (|2Qloh : iLepetit fc Xante ( 2014 ) for experimental observation of BSCs i 


microwave set-ups. Throughout this paper the full destructive interference of two degen¬ 
erate modes resulting in the formation of BSCs in multimodal cavities will be invoked 
to establish two-mode approximation that will be used for finding BSC frequencies and 
shape functions. 


In our recent work Maksimov et al. (12015 ) we showed that the approach known in quan¬ 
tum mechanics as for malism of the effective non-Hermitian Hamiltonian Dittesl ( 2000h : 


Pichugin et al. (|200ll ) could be adapted for solving hard-wall acoustic cavity-duct prob¬ 


lem. The essential feature of the approach is that it allows to recover transmission spectra 
of open-boundary systems relying on the spectral properties of their closed counterparts. 
Technically, the method is the acoustic coupled mode theory in which the pressure field 
within the cavity is represented by the modal variables corresponding to the amplitudes 
of the eigenmodes of th e closed cav ity decoupled from waveguides. We also address the 
reader to the paper by iRacec et all (2009) for a similar approach for scattering in cylin¬ 
drical nanowire heterostructures. 

In this paper, we consider acoustical problem of axisymmetric side-branch cavities in 
two- and three-dimensional ducts of infinite length. Our primary goal is to explore the 
interference phenomena resulting in the formation of BSCs localized within two- and 
three-dimensional ca vities. The BSCs will be analysed in the framework of the acoustic 
coupled mode theory [Maksimov et al. (2 015 ). We will demonstrate that our approach not 
only represents an alternative numerical technique for finding BSCs but also provides 
an efficient analytical tool for construction the BSC shape function in the two-mode 
approximation. In particular, we will see that the two-mode approximation is able to 
a good accuracy reproduce the BSC pressure field within two- and three-dimensional 
acoustic cavities. 


The paper is organized as follows. In Section 2 we briefly overview the acoustic coupled 
mode theory which is our major tool for analysing the scattering problem. The numerical 
techniques for finding BSCs are also discussed. In Section 3 we present our results on 
BSCs in two-dimensional duct-cavity systems. Three-dimensional systems are considered 
in Section 4. Finally, we conclude in Section 5. 






















































































4 A. A. Lyapina, D. N. Maksimov, A. S. Pilipchuk, and A. F. Sadreev, 

2. Acoustic Coupled Mode Theory 

Let us start with the formulation of the problem. Throughout this paper we adopt non- 
dimensionalized quantities. After the time dependance is removed through the substitu¬ 
tion if;(t) = '0exp(— iujt) the pressure field within a cavity-duct structure is controlled 
by the non-dimensional stationary Helmholtz equation 

/ n2 02 \ 

(a? + 5 ? + w)^ + ^ = °- < 21 > 


where uj is the non-dimensional frequency. The x-axis is aligned with the center line 
of the duct-cavity structure. The problem could be reduced to two dimensions under 
assumption that the structure has constant thickness along the z-axis which allows us to 
present the solution of Eq. ED as if = f>(x, y)'ip(z). In what follows we take 'ip(z) = const 
which means that only the first propagating mode is involved. 

Eq. (|2.ip has to be solved subject to Neumann boundary conditions on all sound-hard 
boundaries of the cavity and reflectionless boundary conditions on the cavity-waveguide 
interfaces. There are numerous numerical techniques to impose reflectionless conditions 
on waveguide-cavity interfaces the most efficient of which implement la yers o|_absorbing 
medium with com plex refractivity giving rise to perfectly matched layers \ Berenger ( 19941) 
or complex scaling [Moiseyev ( 19981) methods. For detailed information on t he use of those 
methods for finding acou stic BSCs we address the reader to the papers bv lHein fc Koch 
(120081 ) : \RAn~et ~gj\ (|2Q12h which thoroughly review the subject. In this paper we will us e 
an alternative method, namely, acoustic coupled mode theory [Maksimov et al\ ( 2015 ). 
The method is a recent adaptation of the approach known in quantum mechanics as 
formalism of the effective non-Hermitian Hamiltonian. In this section we only briefly 
review the essentials and introduce the basic concepts that will be needed later on for 
understanding the mechanism of BSC formation. 

To be consistent with the further applications our approach is exemplified on the two- 
dimensional duct-cavity structure depicted in Fig. [l] where we show a two-dimensional 
cavity of length L x and width L y coupled to symmetrically positioned waveguides of 
width d. In what follows we take d = 1. As it was mentioned in the Introduction the 
coupled mode theory relies on the spectral properties of a closed cavity decoupled from 
waveguides. In that context closed means that in the first step the eigenmodes and 
eigenfrequencies of the cavity are computed with the Neumann boundary conditions 


df> 
aa 


= 0 


on 


on the boundary of the resonator dQ including the waveguide-cavity interfaces as well as 
the physical boundaries; a being the local coordinate normal the boundary of the closed 
resonator. In the case of a rectangular cavity the eigenvalue problem is easily solved 
analytically to yield the spectrum as 


uj, 


2 

n,m 


(^) 2 


f tv {m — 1) 

V L y 


2 

; n, m == 1,2,... 


( 2 . 2 ) 


At the same time the eigenmodes are found as products of two factors 




(2 - <*n)( 2 - S m) „„„ (*( n - ViZx + L x) 


Lor.Ln 


cos 


2L, 


cos 


7r(m - 1)(2 y + Ly) 

2 Ln, 


(2.3) 


where 5™, is the Kronecker delta. Notice that the eigenfunctions are unit-normalized, so 
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that 

J ^l,m dxd y = !• (2- 4 ) 

The propagating solutions in the waveguides have a form of plane waves 

< />p ±) (z, y) = ex P (±*fcp^)x P (y), (2-5) 

where symbol +(—) stands for the waves propagating to the right (left), and the functions 
Xp(y) are defined as 

Xp(y) = cos (n —— 1 ^ 2y+ ■ (2.6) 

The dispersion relation reads 

= \A 2 - 71-2 (P - !) 2 ; P = 1,2,.... (2.7) 

The pressure field within the cavity could be constructed as a modal expansion 

V’c = ^ ^ 9n,rn^n,rn • (2.8) 

n,m 


over eigenfunctions Eq. (|2.3|) . The key equation for finding the unknowns g n .m for the 


cavity excited by an incoming wave has the following form Maksimov et al. 


(D-w 2 /)g = -i Yi 

p=l 



VC,pCLC,p 

C=L,R 


(2.9) 
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where / is the identity matrix, the subscript C = L,R specifies the coupling with the left 
and right waveguides, and the unknown vector g contains the coefficients g n ,m- Analo¬ 
gously, the amplitudes ac, p are the coefficients of the modal expansion of the incident 
wave field in the waveguides over the eigenfunctions (|2.5p . while the vector v^p con¬ 
tains the coupling constants between the specific eigenmode (m, n) ([2.3ft and (C,p) 
propagating mode of the waveguides (I231i . Finally, the matrix D has the following form 


D = n 2 


Y ik r Y W G, P ' 

p C=L,R 


( 2 . 10 ) 


where U 2 is a diagonal matrix which carries squared eigenfrequencies u; 2 m Eq. (1231 on 
the main diagonal and the symbol f stands for Hermitian transpose. Notice that the 
coupling to the evanescent modes (p — 1 )tt > uj of the waveguides is accounted for in 
Eq. (j2.10D by the terms with imaginary wavenumbers k v found from Eq. (12.71) . At this 


point we only have to specify the coupling constants According to lPichugin e t al 


(200lj); iMaksimov et al (120151) they are evaluated as overlapping integrals between the 
functions (EU) and the eigenfunctions (1231 on the waveguide-cavity interface 


u C,p 


— J Xp^Pmp 


dr c , 


( 2 . 11 ) 


where rc is the local interface coordinate orthogonal to the wave vector. The generaliza¬ 
tion to the thr ee-d imensional case is straightforward; the coupling constant is a surface 
integral Maksimov et al. (l2015h with rc representing the local coordinate set on a flat 
waveguide-resonator interface. Notice that Eq. (|2.9D for mally conta ins an infinite number 


of unknowns. However, as it was demonstrated in IMaksimov et all (120151 ) for computa¬ 


tional sake the basis can be truncated to a reasonably small number of modes in a given 
frequency window. After the unknown vector g is obtained from Eq. (|2.9j) the coefficients 
of the modal expansion of the outgoing field in the waveguide be, p are found as 

bc, P = ~ac, P + \J 4'7rfc p v ( L p g. (2.12) 

The above described approach is demonstrated to be applicable for fi nding scattering 
functions and reflection amplitudes in two- and three-dimensional cavities IMaksimov et al. 
( 201 5). The goal of the present paper, however, is to adapt it for computing BSCs. The 
BSCs are localized solutions of infinite live-time which exist in the cavity even when no 
incident field is present in the waveguides. Thus, according to Eq. ()2.9|) one obtains an 
eigenvalue problem 

(D - z 2 n I) g n = 0. (2.13) 

Once, an eigenvalue z n = uj n + iq n found from Eq. (|2.13|) is real its eigenvector g n corre¬ 
sponds to an BSC. Unfortunately, Eq. (|2.13D does not represent the standard eigenvalue 
problem because the matrix D Eq. (|2.10|) depends on the spectral parameters k p which 
are linked to the eigenvalue z n through the dispersion relation Eq. (|2.7|) . One of the pos¬ 
sible solutions Sadreev et al\ ( 2006 ) is employing an iterative procedure of solving the 
fix-point equation for real eigenvalues of Eq. (|2.13|) . Alternatively, one of the most effi¬ 
cient techniques for finding complex eigenvalues (|2.13|) is the harmonic inversion method 
which allows to extract the positions and life-times of the resonances from the response 
of the system to an external dri ving. In this paper we do not detail the Harmonic In¬ 
version owing the credit to IWiersig fe Main ( 2008 ) who nicely outlined the method in 
their paper on the fractal Weyl law for chaotic microcavities. In brief, the first step in 
the method is solving Eq. (|2.9|) with some, quite arbitrary, right hand part which, in 
particular, could be a point source placed within the open cavity. The response obtained 
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is then Fourier-transformed to yield a set of nonlinear eq uation for eigenvalues z n which 
is subsequently solved by the use of Pade approximants IWiersig fe Mainl ( 2008 ). 

To conclude this section we retur n to_ th e Friedr i ch-Wi ntgen two-mode interference 
mechanism of BSC formation iFriedrich fe Wintgen (jl985h . Let us check whether it is 
possible that only two modes ifj\ and ^2 of the closed cavity with eigenfrequencies uj 1 
and U 2 contribute to a BSC. Suppose that only the first scattering channel with the wave 
number k\ is open at the frequency of interest. Then, under a further assumptions that 
the coupling to the evanescent channels p > 1 could be neglected we rewrite Eq. (|2.13j) as 
an eigenvalue problem for matrix D with real valued terms uj? arranged along the main 
diagonal in the ascending order 


/ 


D = 


c^o — iki 2 vq 
—iki2voVi 
—iki2voV2 
—iki2voVs 


—iki2voVi 
uj \ — ik\2v\ 
—ik\2v\V2 
—ik\2v\v% 


—ik\2vQV2 
—ik\2v\V2 
uj I — ik\2v\ 
-iki2v 2 vs 


—iki2voVs 
—ik\2v\Vz 
-iki 2 v 2 vz 
ujI — ik\2v\ 


\ 


\ 

, (2-14) 

/ 


where ujq and CJ3 are the eigenfrequencies next to uj\ and UJ 2 , and Vj = v™'™ is the coupling 
constant with the open propagation channel pm 1. Notice, that above we assumed that 
the absolute values of the coupling constants on the left and the right are identical i.e. 
all eigenmodes have compatible symmetry about the y- axis. In fact, the matrix D is split 
into two decoupled subblocks corresponding to the //-symmetric and //-antisymmetric 
eigenmodes so that the //-symmetric and //-antisymmetric problems could be considered 
independently from one another. Under variation of some control parameter, say the 
length of the cavity L x , one can reach the degeneracy point where uq = UJ 2 = ojbsc , 
then the matrix (|2.14[) has a real eigenvalue ujg SC with the eigenvector 


g = 


V 1 V 2 

Tvf 


0, 1 -i/v 2 , 0, ... 


containing only two non-zero entries. Thus, the corresponding BSC shape-function can 
be constructed as 


ipBSC = gU)fpi + 5(2)V’2- (2.15) 

The above equation represents the central mathematical result to be used for analyzing 
BSCs in o pen acoustic resonators. In e ssence, it manifests the well known interference 
mechanism IFriedrich fe Wintgenl (jl985h . In fact, both eigenmodes i/j 1 and ^2 albeit cou¬ 
pled to the waveguides contribute to Eq. (12.15)) in the proportion providing a full de¬ 
structive interference at the degeneracy point so that the resulting BSC function 'ipssc is 
decoupled from the open channel. For details on the realization of the above mechanism 
in the field of quantum mechanics we address the reader to the paper by ISadreev et all 
( 2006 ). We mention in passing that thanks to the symmetry about the duct axis the 
same arguments are valid for the duct-axis antisymmetric modes with &q replaced by k 2 
according to Eq. (12771) . 

The above result Eq. ()2.15j) is exact under assumption that the coupling to the evanes¬ 
cent channels could be neglected. However, we would like to remind the reader that the 
evanescent eigenmodes are present in summation over index p in the right hand part of 
Eq. (|2.10p . Therefore Eq. ()2.15)) is not an exact solution, and could only be seen as the 
basis for a two-mode approximation for BSCs occurring at the degeneracy points. Some 
remarks are due on the validity of the two-mode approximation. As it will be seen later on 
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P \(n,m) 

( 1 , 1 ) 

(1,2) 

(1,3) 

(1,4) 

(4,1) 

(4,2) 

(2,3) 

(2,4) 

1 

0.408 

0 

-0.367 

0 

0.577 

0 

-0.520 

0 

2 

0 

0.245 

0 

-0.441 

0 

0.346 

0 

-0.624 

3 

0 

0 

0.173 

0 

0 

0 

0.245 

0 

4 

0 

0.021 

0 

0.082 

0 

0.030 

0 

0.116 

5 

0 

0 

0.034 

0 

0 

0 

0.049 

0 

6 

0 

0.007 

0 

0.024 

0 

0.011 

0 

0.034 


Table 1. Numerical values of coupling constants v2’™ at L x = 3, L y = 2. 


due to the structure of Eq. (j2.10D the evanescent coupling results in the presence of other 
than the two basic modes in the modal expansion of the BSC shape function, as well as 
in a shift of the true BSC point from the degeneracy point. Nevertheless, one can argue 
that the coupling to the evanescent channels is small relative to the coupling to the open 
channel as we consider the eigenmodes residing in the frequency window where only one 
symmetric and one antisymmetric propagating channels are open. In order to illustrate 
this we collected in Table [T] the numerical values of the coupling consta nts for so me rel¬ 
evant eigenmodes evaluated on the basis of the analytical results from iMaksimov et al. 
(1201SL One can see from Table [T| that the coupling constants rapidly decay with the 
growth of the channel number p. In fact, in can be demonstrated that any coupling 
constant decays away from its maximal value with the asymptotic law ~ 1/p 2 . 

Finally, notice that the two-mode approximation has no restrictions related to the 
size of the cavity. One can see from Eq. (|2.14|) that once a degeneracy occurs in the 
closed cavity it always results in the formation of a BSC through the Friedrich-Wintgen 
mechanism irrelevant to the values of L x , and L y . If the size of the cavity is isometrically 
increased whilst the width of the waveguide is kept unchanged, the next to the BSC 
eigenfrequencies would shift closer to the BSC eigenfrequency. The real valued 

diagonal term in Eq. (|2.13[) will drop as (cJq 3 ~ ^b sc) ~ 1/5? S' = L x L y because the 
mean distance between the eigenfrequencies is inversely proportional to the width/length 
of the cavity. On the other hand the coupling constants Vj in matrix D (|2.14j) will drop as 
Vj ~ 1/VS because the eigenfunctions Eq. (|2.3[) contributing to the overlapping integrals 
Eq. (|2.11|) have the normalization constant y/S in denominator. In consequence of this 
the off-diagonal terms VjVj> in matrix ([2.14)) drop as VjVj> ~ 1/S. Therefore, although 
in larger cavities the next to the BSC eigenvalues appear closer to the BSC point their 
coupling to the basic BSC modes through the evanescent channels drops at the same 
rate as the resonant term cjq 3 — ^bsc- This leads us to a conclusion that the two-mode 
approximation does non break down as the size of the cavity is increased. Notice, that 
the same holds true in the three-dimensional case where the coupling constant drops as 
Vj ^ 1/y/v with V as the volume of the cavity. 

3. Bound states in two dimensional duct-cavity structure 

We start with computing the transmittance T = |^i 5 J r | 2 for the two-dimensional cavity 
in Fig. [T| under assumption that an incident wave in the first propagation band p = 1 
enters the cavity from the left waveguide. For solving the scattering problem we used the 
approach described in Section 2 . The numerical vdues_ of_the _coupling_constants v//™ 
in Eq. (j2.9l2.10p are taken from our previous work IMaksimov et all (j2015t h In Fig. [2] we 
demonstrate the dependence of the transmittance T on the frequency of the incoming 
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1.0 

0.8 

0.6 

0.4 

0.2 

0.0 

Figure 2. (Color on-line) Transmittance for a two-dimensional cavity of length L y — 2. Only 
the first symmetric channel with p — 1 is open. The dashed lines show eigenfrequencies of the 
modes symmetric about the y- axis; solid line - antisymmetric ones. The positions of the BSCs are 
shown by open circles and stars for y-symmetric and ^-antisymmetric modes, correspondingly. 



wave uo and the length of the cavity L x . Notice that the frequency varies in the range 
[0, 27r] which means according to Eq. (|2.7|) that only one symmetric scattering channel is 
open in each waveguide. In Fig. [2] we also depict the frequencies of the eigenmodes of the 
closed cavity symmetric about the duct axis Eq. (|2.2j) m = 1,3, 5, — The eigenfrequen¬ 
cies which correspond to the modes symmetric about the y- axis, later to be referred to 
as ?/-symmetric, are shown by dashed lines. The solid lines show the eigenfrequencies of 
?/-antisymmetric eigenmodes. One can see that the degeneracy points of eigenfrequencies 
with the same symmetry about the y -axis adjoin to an abrupt change of the transmit¬ 
tance; the transmittance drops from unit to zero as the frequency is swept across the 
point where the collapse of Fano resonance occurs. The complex resonance frequencies 
for the cavity are computed by harmonic inversion method. In the given range of pa¬ 
rameters we found five BSCs. Three of them originate from the modes symmetric about 
the 7/-axis (shown by circles) and two from the antisymmetric ones (shown by stars). 
For further convenience we numerate the BSCs from 1 to 5. The complex eigenfrequen¬ 
cies are plotted against the length of the resonator L x in Fig. [3] where one can see that 
the imaginary part y n turns to zero as one approaches the BSC point. In Fig. 0] we 
demonstrate the numerical solution for BSC 3 along with the corresponding coefficients 
of the modal expansion (Oil . One can clearly see that only two eigenmodes dominate 
in the BSC function 'ipssc • The ^-antisymmetric BSC 3 in Fig. 0] is formed by eigen¬ 
modes Vh = Vt,i { x ,y), and ^2 = V^ 2 , 3 y) with the corresponding coupling constants 
v\ = v^ \ = —Vr\ = 0.600, and V 2 = v^\ = —v^\ = —0.540. To assess the accuracy of 
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L L 

X X 

Figure 3. (Color on-line) Complex eigenfrequencies against the length of the two-dimensional 
resonator L x \ real parts - left panel; imaginary parts - right panel. The //-symmetric resonances 
are shown by dash brown lines; //-antisymmetric by blue solid lines. The positions of BSCs from 
Fig.m are shown by red open circles and stars for //-symmetric and //-antisymmetric modes, 
correspondingly. 

the two-mode approximation we introduce parameters 


Aw= 

(3.1) 

UBSC 

1 - t* 1 

A l x = 1 - ai , 

-L'x 

(3.2) 


with u/*, and L* as the frequency and the length of the cavity corresponding to the 
degeneracy point in which the BSC would occur according the two-mode approximation 
Eq. (|2.14|) . For BSC 3 we obtained A uj = 0.019, and A L x = 0.018. In Fig. [4] we also 
plot the BSC shape function constructed by Eq. (|2.15|) . One can see from Fig. 0 that the 
two-mode approximation reproduces our numerical result to a good accuracy. We also 
found that the same holds true for the other four BSCs in Fig. [2j Therefore, the other 
BSC functions are not presented in this paper. 

The same procedure was repeated in case of the incoming wave in the first antisymmet¬ 
ric channel for cavities with L y = 2. As before we observed that the two-mode approxima¬ 
tion is valid in the situation when only one antisymmetric channel is open in each waveg¬ 
uide. For brevity we omit a detailed description of our findings restricting ourselves to 
one example of a duct axis antisymmetric BSC shown in Fig. [5] The //-symmetric BSC in 
Fig. [5] is formed by eigenmodes t/q = ^ 5 , 2 OC //), and ^2 = ^p 3 , 4 ( 2 ?, y) with the correspond¬ 
ing coupling constants v\ = = v^\ = —0.692, and V 2 = v^\ = v^\ = 0.384. For 

the corresponding relative errors Eqs. (|3.1I3.2[) we found A uo = 0.004, and A L x = 0.004. 
As it was mentioned in Section 2, we expect that the two-mode approximation could 
be applied irrelevant to the cavity size. This is confirmed by numerical tests for larger 
cavities L y = 4. As an example in Fig. [6] we present the duct-axis antisymmetric BSC 
formed by eigenmodes ^1 = ^5,8 (A2z)> and ^2 = ^ 7,6 y) with the coupling constants 
v\ = = —0.371, and V 2 = v 7 £\ = v 7 jf\ = 0.402. The values of the relative 
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Figure 4. (Color on-line) BSC Asymmetric about the duct axis of a two-dimensional duct-cavity 
structure uobsc — 3.396, L y = 2, L x = 2.776; a) Pressure field of the numerically computed BSC 
shape function; b) Absolute values of the modal coefficient of the numerically computed pressure 
field; c) BSC shape function found from the two-mode approximation Eq. (I2.15D . 



Figure 5. (Color on-line) Bound state in the continuum antisymmetric about the duct axis of 
a two-dimensional duct-cavity structure ujbsc = 5.385, L y — 2, L x = 2.441; a) Pressure field 
of the numerically computed BSC shape function; b) Absolute values of the modal coefficient 
of the numerically computed pressure field; c) BSC shape function found from the two-mode 
approximation Eq. iimn . 
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Figure 6 . (Color on-line) Bound state in the continuum antisymmetric about the duct axis of 
a two-dimensional duct-cavity structure ujbsc = 6.502, L y — 4, L x = 3.627; a) Pressure field 
of the numerically computed BSC shape function; b) Absolute values of the modal coefficient 
of the numerically computed pressure field; c) BSC shape function found from the two-mode 
approximation Eq. (I2.15D . 


errors A u = 0.003, and A L x = 0.007 indicate that the two-mode approximation holds 
true for larger cavities. 

As it was mentioned in the introduction the occurrence of BSC is associated with a 
collapse of Fano resonance. To illustrate this in Fig. [7| we plot the dependance trans¬ 
mittance vs. frequency for three different len gth s of th e reso nator L x close to the BSCs 
point. It was recently demonstrated bv lBulgakov &; Sadreev ( 20141 ) that if the BSC point 
is approached in the parametric space along the certain path corresponding to the peak 
of the Fano resonance, the long-lived complex mode related to the BSC provides the 
dominant contribution to the scattering function. In consequence of this the scattering 
function near the BSC point bears the maximal resemblance of the BSC if the parameters 
are tuned to the peak of the Fano resonance. This is illustrated in Fig. [8] where one can 
see that the real part of the scattering function is almost identical to the BSC from Fig. 
|4] meanwhile the imaginary is present as a background providing the net power flux in 
the system. Thus, one can experimentally observe a signature of BSC in the form of a 
large amplitude acoustic resonance even though in a realistic set-up the parameters are 
always detuned from the true BSC point due, for example, to fabrication inaccuracies. 


4. Bound states in the continuum in three-dimensional open cavities 

Let us consider the axisymmetric cylindrical resonator with two coaxially attached 
cylindrical waveguides of radius a = 1 as depicted in Fig. [9ja). After the azimuthal 
variable is separated we can restrict ourself to the case of zero angular momentum bearing 
in mind that our results can be easily generalized to the non-zero case taking into account 
that by virtue of Eq. (|2.11|) only the modes of the same rotational number are coupled 
on the waveguide-cavity interface. In our case the cavity eigenfrequencies are given by 
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Figure 7. (Color on-line) Collapse of Fano resonance in the vicinity of a BSC point under 
variation of the length of the resonator L x . a) BSC 4 in Fig .\4\ujbsc — 3.396, L y = 2, L x = 2.776. 
Resonator lengths: L x = 2.75 - solid blue line, L x = 2.70 - dashed green line, L x = 2.60 - dash-dot 
red line, b) Duct axis antisymmetric BSC in Fig.^cassc = 3.384, L y — 2,L X = 2.441. Resonator 
lengths: L x = 2.460 - solid blue line, L x = 2.485 - dashed green line, L x = 2.510 - dash-dot red 
line. 
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Figure 8. (Color on-line) Scattering function at the peak uj — 3.410 of the Fano resonance Fig. 
[7| related to BSC 3 shown in Fig. [H L y = 2, L x = 2.75. a) Real part of the scattering function, 
b) Imaginary part of the scattering function. 
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Figure 9. Bound states in axisymmetric duct-cavity structure, a) Cylindrical resonator of radius 
R and length L x with two coaxially attached waveguides of radius a — 1. b) Pressure field 
isosurfaces of BSC 2 in Fig. fl2l symmetric with respect to the central section of the resonator; 
dark blue 'ipBSC — —0.1, light red 'ipBSC — 0.1. c) Pressure field isosurfaces of BSC 7 in Fig. 
m antisymmetric with respect to the central section of the resonator; dark blue ^bsc = —0.1, 
light red 'ipBSC = 0.1. 


uj, 


2 

m,n 


\f4n , ^(n-l) 2 ] 

.R 2 L l 


(4.1) 


where R and L x are the radius and length of the cylindrical resonator, and is the 
m-th root of the equation Jo(/im) = Ji(Rm) = 0,m = 1,2 ,... for the derivative of the 
first order Bessel function. The corresponding eigenfunctions are 


4>m,n{r,x) 


1 

^nRJ 0 {n m ) 




(4.2) 
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Figure 10. (Color on-line) Transmittance for a cylindrical resonator with the radius R = 2 with 
two coaxially attached cylindrical waveguides of radius a — 1 vs. frequency uj and the length 
of the resonator L x . The dashed lines show eigenfrequencies of the modes symmetric about the 
yOz-plane; solid line - antisymmetric ones. The positions of the BSCs are shown by open circles 
and stars for ^/Oz-symmetric and yOz-antisymmetric modes, correspondingly. 


where r is the radial coordinate in the T/Oz-plane while the axial modes are given by 


l/>n(x) = 


SI 


■ COS 


7T (n — l)(2x + L x ) 
2L X 


, n = 1,2 ,... 


(4.3) 


As before th e transm ission problem was solved by the use of the acoustic coupled mode 
theory [Maksimov et all ( 20151) and the complex eigenvalues were obtained by the har¬ 
monic inversion method. The frequency uj varies in the range [0,/i2],/i2 = 3.831 which 
means that only one propagating channel is open in each waveguide. The transmittance 
along with the spectrum of the closed resonator and the positions of BSCs is shown in Fig. 
m The complex eigenfrequencies are plotted in Fig. [lT] against the length of the resonator 
L x . In Figs. II21451 we demonstrate the shape functions of two BSCs occurring through the 
two mode interference mechanism. BSC 2 in Fig.[12]with A uo = 0.012, and A L x = 0.012 
is formed by eigenmodes Vh = and ^2 = ^ 3,1 (r, x) with the corresponding 

coupling constants v\ = v£ 1 = 1 = 0.310, and V 2 “ v£ 1 = 1 = 0.289, while yOz- 

antisymmetric BSC 7depicted in Fig.[l5]is formed by ?/q = ^ 2 , 4 (r, x), and ^2 = ^ 4 , 3 (t*, x) 
with the coupling constants v\ = vp x = —0.078, and V 2 = 1 = ~ v r 1 = 0.389. 

A three dimensional visualization of BSCs from Figs. 1121131 is offered in Fig. [9] b,c) in 
form of pressure field isosurfaces. One can see that of all BSCs in Fig. [10] BSC 7 deviates 
the most from the corresponding degeneracy point A uj = 0.016, and A L x = 0.066. The 
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Figure 11. (Color on-line) Complex eigenfrequencies against the length of the three-dimensional 
resonator L x \ real parts - left panel, imaginary parts - right panel. The yOz-symmetric resonances 
are shown by dash brown lines; yOz-antisymmetric by blue solid lines. The positions of BSCs 
from Fig. [10] are shown by red open circles and stars for yOz-symmetric and yOz-antisymmetric 
modes, correspondingly. 


reason for this can be understood from Figs. l9lc). [T3l where one can clearly see that BSC 
7 extends far into the waveguides due to the coupling to the evanescent modes which is 
neglected by the two mode approximation. One can conclude that the coupling to the 
evanescent modes becomes important just below the next channel propagation thresh¬ 
old as the channel function extends far into the waveguide. Nevertheless, the two-mode 
approximation is still valid in the three-dimensional case reproducing the numerically 
computed shape functions to a good accuracy. 

In this paper, we restricted ourself to coaxial systems in which the azimuthal num¬ 
ber is a constant of motion so that the rotational vibration modes do not participate 
in formation of BSC depicted in Figs. [12] [13] It is worth noticing, however, that the 
cavity eigenmodes with azimuthal number v ^ 0 are symmetry protected BSCs as their 
eigenfrequencies are below the u-th propagating threshold in cylindrical waveguides. 


5. Conclusion 

We considered sound transmission through two- and three-dimensional acoustic res¬ 
onators and demonstrated the existence of singular points in which the unit transmittance 
coalesces with the zero transmittance in the points of collapse of Fano resonances. These 
singular points occur due to a full destructive interference of two degenerate modes of 
the same symmetry. As a result a certain superposition of these two modes is a trapped 
solution corresponding to a bound state in the continuum or embedded trapped mode 
localized within t he resonat or. The above mechanism of wave localization was first pro¬ 
posed bv lFriedrich fe Wintgenl (1 19851). and has been so far e xperi mentally realized only 
in microwave set-ups Lepetit et al | ([201 Oh : iLepetit &; Kant^l (2014). In the present paper 
we invited the reader’s attention to acoustic cylindrical resonators coaxially connected to 
semi-infinite waveguides where the degeneracy points could be accessed in the parametric 
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Figure 12. (Color on-line) BSC 2 in a three-dimensional duct-cavity structure 
ujbsc — 2.367, R — 2,L X = 2.655; a) xOy-plane section of the pressure field of the numerically 
computed BSC shape fnction; b) Absolute values of the modal coefficient of the numerically 
computed pressure held; c) BSC shape function found from the two-mode approximation Eq. 
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Figure 13. (Color on-line) BSC 7 in a three-dimensional duct-cavity structure 
ujbsc — 3.557, R — 2,L X = 3.387; a) x(h/-plane section of the pressure held of the numeri¬ 
cally BSC; b) Absolute values of the modal coefficient of the numerically computed pressure 
held; c) BSC shape function found from the two-mode approximation Eq. jnsj. 
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space under variation of the length of the resonator. We speculate that such a tuning 
could be performed in a realistic acoustic experiment by the use of piston-like hollow-stem 
waveguides tightly fit to the interior boundaries of a cylindric cavity. For analysing the 
boun d states in t he continuum we successfully employed the acoustic coupled mode the¬ 
ory iMakshno^^L^ ( 2015 1 which allowed us to analytically predi ct the eigenfrequen cies 
of the bound states as well as their shape functions. Until recently Hein et al. (2012) the 
bound states in the continuum were independently studied in physics and fluid mechanics 
with little to none cross-fertilization. In the present paper we attempted to bridge that 
gap to the mutual benefit of the both communities. 

Acknowledgement: This work has been supported by Russian Scientific Foundation 
through Grant 14-12-00266. 
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